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Abstract 

The one-dimensional effective-mass Klein-Gordon equation for the real, and non-PT- 
symmetric/non-Hermitian generalized Morse potential is solved by taking a series expansion for 
the wave function. The energy eigenvalues, and the corresponding eigenfunctions are obtained. 
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I. INTRODUCTION 



In the past few decades, there has been an increased interest to find the exact solutions 
of the non-relativistic and relativistic equations. The solutions of the Schrodinger equation 
has been studied by using different methods based on perturbative and non-perturbative 
approaches [1-5]. The Klein-Gordon (KG) and Dirac equations have been also studied for 
different type of potentials such as Aharonov-Bohm (AB) potential [6], the AB plus Dirac 
monopole potential [7, 1 1] , generalized Hulthen, harmonic, and linear potentials, generalized 
asymmetrical Hartmann potentials, for a uniform magnetic field, pseudoharmonic oscillator, 
and exponential-type potentials [12-20]. 

Recently, there has been considerable work on the PT-symmetric quantum mechanics. 
Following the initially studies of Bender et al. [21], the PT-symmetric theory has been 
successfully studied by many authors because of their importance, and also for using in 
different research areas, such as quantum field theory, and nuclear physics [22, 23]. In the 
present study, we take into account the form of the generalized Morse potential having 
also non-PT symmetric property to achieve the solutions within the posit ion- dependent 
mass formalism, and study the effects of the mass varying with spatially coordinate on the 
solutions of the case of constant mass. It could be interested to solve the relativistic KG 
equation in this point of view. 

The solutions of the above wave equations are taken in general for constant or time- 
dependent masses [24, 25]. There are also many examples of physically important systems, 
for which the mass depends on coordinates. So far, solving the above equations for the 
case of spatially dependent mass has been studied by many authors [26-34] for different 
types of the mass distributions such as an exponentially, and quadratic mass functions [27], 
inverse-quadratic dependence of mass [35], trigonometric mass-distributions [36], and mass 
functions of the forms m(r)=r a , especially for three-dimensional problems [36, 37]. We will 
use also a mass distribution having an exponential form, and study in the half-plane to 
achieve a physically acceptable (finite) mass function, which has a decreasingly behavior in 
this range. 

In the present work, we study the bound state solutions of the KG equation with real 
and non-PT symmetric/non-Hermitian forms of generalized Morse potential in the case of 
a co ordinate- dependence mass distribution. In order to find the energy spectra, and the 
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corresponding wave functions we prefer to use an infinite series for the wave function. This 
approach is a powerful technique to solve second order differential equations especially for the 
potential forms having two and/or more terms [38, 39]. We obtain also the energy spectrum, 
and the corresponding eigenfunctions in the case of constant mass for two potential types. 

The organization of the work as follows. In section II, we obtain the exact bound-state 
energy spectra for real, and non-PT symmetric/non-Hermitian cases of generalized Morse 
potential, and the corresponding eigenfunctions by taking an infinite series for the wave 
function in the case of a convenient mass distribution. We summarize the numerical results 
in one table, and two figures to see the effects of the mass depending on coordinate. For this 
aim we do the numerical analysis for three different diatomic molecules such as H 2 , LiH, and 
HC1. We give also the results for the case of constant mass, and present our conclusions in 
Section III. 



II. BOUND STATE SOLUTIONS 
A. Real Generalized Morse Case 



The one-dimensional KG equation is written in terms of scalar and vector potential as 
following [40] 



d 2 i/j(x) 1 



[E - V v (x)f - [mc 2 + V s {x)} 2 i){x) =0, (1) 



dx 2 h 2 c 2 

where m is the mass of the particle, E denotes the energy, and V s (x) and V v (x) are the 
scalar and vector parts of the potential, respectively, and c is the velocity of light. In the 
absence of the vector potential, the KG equation can be written as 



d 2 ^(x) 2m ( V 2 (x) , >. i 



da; 2 fi 2 \ 2mc 2 
which is a Schrodinger-like equation with the energy eigenvalue e = — m ^-. 

We consider the scalar potential V s (x) in Eq. (2) as generalized Morse potential, which 
can be used to describe the vibrations of a two-atomic molecule, as the following 



V s (x) = V^- 2 ^ - V 2 e~P x , (0 < x < +oo) 



(3) 



where Vi, and V2 are two real parameters, and x = (r — r )/r . The parameter (3 is ar , 
here, a is the potential width, and r is the equilibrium distance. 
Here, we prefer to use the following mass-distribution 



m(x) = m + rriie ^ x , (4) 

where m , and mi are real parameters. This distribution is finite at infinity, and enables us 
to analyze the results in the case of constant mass. 

Substituting Eqs. (3) and (4) into Eq. (2), than we have 



+ {j^T2^ 2 - m oC 4 ] + ^ [V 2 m c 2 - m m lC 4 ]e 

- J-7 Wi + 2Vim c 2 - 2V 2 m lC 2 + m\c A \^ 2fix 
n c 2 

+ Jr~ 2 [ViV 2 - V imi c 2 }e-^ - J£ e-V*U{x) = , (5) 

ft C ft c * 

By using the transformations on coordinate and wave function 



z = e- fix , #s)=-^=0(z), (6) 
and with the help of the following parameters 



A 1 = Q 2 (E 2 -m 2 c 4 ) + i, 

A 2 = 2Q 2 (V 2 m c 2 - m mic 4 ) , 

^3 = Q\V 2 + 2V imo c 2 - 2V 2 m lC 2 + mlc 4 ) , 

A 4 = 2Q 2 (V 1 V 2 - V imi c 2 ) , 

A 5 = Q 2 V 2 . (7) 

we get 

^ + (~A 3 + ^ + ^ - A,z- A 5 A <p(z) = 0. (8) 
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where Q 2 = l/h 2 c 2 /3 2 in the above expressions (0 < z < 1). It would be interested the 
results obtained from the last equation for a special case A 2 =A 4 = 0, which gives 



^ + - ^ 2 ) *{z) = A 3 <j>(z) , (9) 

This equation has a similar form with the ones given by Eq. (7) in Ref. [17] for A 5 — > a 2 , 
Ai ->■ —£'(£'+ 1), and A 3 ->■ -A (a,£', and A are the parameters used in Ref. [17]). Eq. (9) 
corresponds to the equation of the harmonic oscillator with centrifugal potential barrier, so 
we could give the solutions as 



r 4r 2 P\ o 1 -i 1/2 

e n = T{mlc 4 + -^(mic 2 + m ) 2 [2n + 1 + ^ - - IfQ 2 ] - ^} . (10) 

We turn to the solution of Eq. (8), and write the wave function as [38, 39] 



c/,( z ) = tf'+W)** J2 a n z 2n+L+1 ' 2 . (11) 

n=0 

Substituting Eq. (11) into Eq. (8), and equating of the coefficients to zero, we get the 
following identities among coefficients, and the expressions for p, and q 



X n a n + Y n+ ia n+ i + Z n+ 2d n+ 2 — , (12) 

where 



X n — 



Y„. = 



2q{2n + L + 1) - A 3 , 
A 2 + p(4n + 2L + l), 
4n(n + L) + 2L 2 . 



(13) 
(14) 
(15) 



and 



q 2 = A 5 
2pq = A 4 . 
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(16) 
(17) 



where the new parameter L is defined as L 2 — > Q 2 (E 2 — m^c 4 ) , and we choose the parameters 
as p= Qiy-i — mic 2 ), and q= —QV\ to obtain a physically solution. 

On the other hand, X„, Y n and Z n must satisfy the following condition to determine the 
coefficients in the system of equations given in Eqs. (13)-(15) 



det 



Y Z 1 ... 

Xo Yi Z 2 
X± Y 2 Z 3 









X n _x Y n 



= 0, 



(18) 



In the case of a n ^ 0, but a n+ i =a. n +2 
leads to the following algebraic equation 



in Eq. (12), we impose X„ = 0. This 



2q{2n + L + 1) = A 3 . 



(19) 



Substituting the values of the parameters defined in Eq. (7) into Eq. (19), and by using 
p, and q obtained from Eqs. (16)-(17), we have the energy eigenvalues of the generalized 
Morse potential as 



E n = ±{ m 2 c 4 + 



where 



Q 2 



2n + 1 + QV 1 



+ 



mic 2 V2 



1/2 



2(2n + 1) + 2QV 1 + Qm lC 2 V 2 



(20) 



V 2 



Vo 



■2 , • "2 . rriiC 



(21) 



It is seen that the energy levels of the particle and antiparticles are split around zero. 
On the other hand, the energy spectrum of the generalized Morse potential for the position- 
dependent mass is real. The part of the energy eigenvalues coming from the coordinate 
dependence of the mass is dependent on the quantum number n, and also on the potential 
parameters Vi, and V 2 . 
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The corresponding wave functions can be written as the following 



4> n (z) = {a z 1/2 + a lZ 5/2 + . . )z L exp Q(V 2 - m lC 2 )z 



QV l 



(22) 



where a;(i = 0, 1 . . .) are the coefficients of the series. 

In order to obtain the energy spectrum of the real generalized Morse potential for the 
case of constant mass, we have to set mi = in Eq. (20), than we get 



According to Eq. (23), the energy levels of the particle and antiparticles are symmetric 
about the zero in the limit of the constant mass. On the other hand, the energy spectrum 
of the real generalized Morse potential is purely real as in the first case. 

We summarize the numerical results in Table I for the H 2 , LiH, and HC1 molecules in 
the case of the constant mass m obtained from Eq. (23). We give the bound state energies 
of the particles (+E n ), and antiparticles (-E n ) for different n- values in Table I. Further, we 
plot the dependence of the ground states of the above molecules on mi for H 2 , and LiH, 
HC1 molecules in Fig.l, and Fig.2, respectively. We do the numerical analysis in the range 
10~ 6 < < 10~ 4 , where M=m /mi, and (p) denotes the particle case, and (a) denotes the 
antiparticle one in figures. 

The corresponding eigenfunctions in the case of the constant mass can be obtained from 
Eq. (22) by setting m 1 = 0. 

B. Non-PT Symmetric and non-Hermitian Generalized Morse Case 

Now we consider the complex case of the generalized Morse potential, where the potential 
parameters are written as 




(23) 



Vi = v 2 — v% + 2iviv 2 , 

V 2 = v± + iv 2 + 2v 3 (vi + iv 2 ) . 



(24) 
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where v[s(i = 1,2,3) are the real parameters, and (3 = 1. This complex form of the 
generalized Morse potential is called non-PT symmetric/no n-Hermitian potential form of 
the potential [41, 42]. 

Substituting Eq. (24) into Eq. (21), than we have from Eq. (20) for the energy spectra 
of the non-PT symmetric/non-Hermitian generalized Morse potential the following 



E„ 



± 



2„4 



m c 



+ 



Q 



12 



2n + 1 + Q'V{ 



+ 



1/2 



2(2n + 1) + 2Q'V{ + Q'm lC 2 V 2 ' 



(25) 



where 



~ 2(v 3 + 1) 2 2 - 2v 3 + l m x c 2 1 

2 vi + w 2 2(vi + w 2 y he 2 

The corresponding eigenfunctions can be written as 

<j> n (z) = (a Q z L+l ' 2 + a lZ L+ ^ 2 + ...) exp 
The energy spectrum for the case of constant mass is given as 

£r=° = ± fan + 1 + Q'VA \ + { ^ - V, (28) 

V V \/ Q' 2 \2n + l + Q'Vj 

It is clearly seen that the energy eigenvalues of the non-PT symmetric/non-Hermitian 
generalized Morse potential has real and imaginary parts, which is coming from the coor- 
dinate dependence of the mass, but in the constant mass limit the energy spectrum is real. 
On the other hand, the energy levels of the non-PT symmetric/non-Hermitian generalized 
Morse potential is symmetric about zero as in the real potential case. The eigenfunctions 
for the case of constant mass are obtained by setting mi = in Eq. (27). 



Q(V 2 - m lC 2 )z 



QVi 2 

— — z 



(27) 



III. CONCLUSION 



We have obtained the solutions of the KG equation for the real and non-PT 
symmetric/non-Hermitian generalized Morse potential with position dependent mass. The 
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energy spectra and the corresponding wave functions have been obtained by using an infi- 
nite series for the wave function. We have found that the real generalized Morse potential 
has a real energy eigenvalues in the case of the position-dependent mass. We have sum- 
marized first the energy eigenvalues of three diatomic molecules for different values of n, 
and than plot two different graphs where the energy eigenvalues varying with mass to see 
the effects of the mass depending on spatially coordinate. We have also studied the non-PT 
symmetric/non-Hermitian case of the potential, and pointed out that the energy eigenvalues 
are imaginary in this case. We have found that the energy levels are purely real in the case 
of the constant mass as given in Eqs. (23), and (28) for both form of the generalized Morse 
potential. 



IV. ACKNOWLEDGMENTS 



This research was partially supported by the Scientific and Technical Research Council 
of Turkey. 



9 



J. P. Killingbeck, A. Grosjean, and G. Jolicard, J. Phys. A: Math. Gen. 34, 8309 (2001). 
A. Arda, Turk. J. Phys. 28, 223 (2004). 

J. Chen, L. C. Kweck, and C. H. Oh, Phys. Rev. A 67, 012101 (2003). 
S. N. Biswas et al., J. Math. Phys. 14, 1190 (1973). 

A. Khare, and U. P. Sukhatme, J. Phys. A 26, L901-L904 (1993) [arXiv: hep-th/9212147| . 



C. R. Hagen, Phys. Rev. D 48, 5935 (1993), [arXiv: [hep-th/9308009| 
V. M. Villalba, [arXiv: [hepjth/9409102 1. 



A. S.de Castro, Int. J. Mod. Phys. A 22, 2609 (2007) [arXiv: [hep^th/0511010| 



A. S.de Castro, Phys. Lett. A 338, 81 (2005), [arXiv: (hepjth/0502201 1 



A. S.de Castro, Phys. Lett. A 346, 71 (2005), [arXiv: (hepjth/05072181. 



A. S.de Castro, Phys. Lett. A 342, 53 (2005) [arXiv: hep-th/0409216] . 
W. C. Qiang, R.-S. Zhou, and Y. Gao, Phys. Lett. A 371, 201 (2007). 
N. Saad, [arXiv: math-ph/0709.4014]. 

R. Giachetti, and E. Sorace, Phys. Rev. Lett. 101, 190401 (2008) [arXiv: hep-th/0706.0127]. 

K. Bhattacharya, [arXiv: hep-th/0705.4275]. 

A. S. Dutra, and M. B. Hott, [arXiv: quant-ph/0705.3447]. 

C. Gang, C. Zi-Dong, and L. Zhi-Mei, Chi. Phys. 13, 279 (2004). 

C. Gang, Phys. Lett. A 339, 300 (2005). 

C. Gang, Acta. Phys. Sin. 50, 1651 (2001). 

L. Z. Yi, Y. F. Diao, J. Y. Liu, and C. S. Jia, Phys. Lett. A 333, 212 (2004). 

C. M . Bender, S. Boettcher, and P. N. Meisinger, J. Math. Phys. 40, 2201 (1999). 

S. M. Ikhdair, and R. Sever, Ann. Phys. 16, 218 (2007). 

C. S. Jia, P. Q. Wang, J. Y. Liu, and S. He, Int. J. Theor. Phys. 47, 2513 (2008), and references 
therein. 

A. Mostafazadeh, J. Phys. A 31, 6495 (1998). 
A. Mostafazadeh, Phys. Rev. A 55, 4084 (1997). 
A. D. Alhaidari, Phys. Lett. A 322, 72 (2004). 

A. S. Dutra, and C. A. S. Almeida, Phys. Lett. A 275, 25 (2000). 

B. Gonul, B. Gonul, D. Tutcu, and O. Ozer, Mod. Phys. Lett. A 17, 2057 (2002). 



10 



B. Gonul, O. Ozer, B. Gonul, and F. Uzgun, Mod. Phys. Lett. A 17, 2453 (2002). 
B. Gonul, and M. Kocak, Chin. Phys. Lett. 20, 2742 (2005). 



B. Gonul, and M. Kocak, [arXiv: quant-ph/05 12035 1. 



C. Tezcan, and R. Sever, J. Math. Chem. 42, 387 (2007) [arXiv: |quant-p h/0604041|. 

I. O. Vakarchuk, J. Phys. A: Math. Gen. 38, 4727 (2005). 

C. Gang, and Z. Chen, Phys. Lett. A 331, 312 (2004). 

L. Jiang, L.-Z. Yi, and C.-S. Jua, Phys. Lett. A 345, 249 (2005). 

A. D. Alhaidari, Phys. Rev. A 66, 042116 (2002). 

G. X. Ju, Y. Xiang, and Z. Z. Ren, [arXiv: |quant-ph/0601005] . 



S. H. Dong, Z. Ma, and G. Espozito, Found. Phys. Lett. 12, 465 (1999). 

S. H. Dong, Int. J. Theor. Phys. 39, 1119 (2000); 40, 569 (2001). 

A. S. de Castro, Phys. Lett. A 338, 81 (2005). 

Z. Ahmed, Phys. Lett. A 282, 343 (2001). 

G. Bagchi, and C. Quesne, Phys. Lett. A 273, 285 (2000). 

E. D. Filho, and R. M. Ricotta, Phys. Lett. A 269, 269 (2000) [arXiv: |hep-th/9910254] . 



11 



TABLE I: The dependence of bound states on n in MeV for H 2 (D=38267.76 cm -1 , a=l,9426l _1 , 
r = 0,74161, and m = 0, 50391amu), LiH (D=20287 cm" 1 , a=l,1280A _1 , r = 1,59561, and 
m = 0,8801221 amu), and HC1 molecule (D=37255 cm -1 , a=l,8677l _1 , r = 1,27461, and 
m = 0, 9801045 amu) [43] with £ = 0. 
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FIG. 1: The variation of ground states of the H2 molecule with mi. 
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FIG. 2: The variation of ground states of the LiH, and HCI molecules with mi. 
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